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Abstract. Mass loss is a determinant factor which strongly affects the evolution and the fate of 

massive stars. At low metallicity, stars are supposed to rotate faster than at the solar one. This 

favors the existence of stars near the critical velocity. In this rotation regime, the deformation of 

the stellar surface becomes important, and wind anisotropy develops. Polar winds are expected 

to be dominant for fast rotating hot stars. 

QQ " These polar winds allow the star to lose large quantities of mass and still retain a high 

f^ . angular momentum, and they modifie the evolution of the surface velocity and the final angular 

f*^ ' momentum kept in the star's core. We show here how these winds affect the final stages of 

fSJ , massive stars, according to our knowledge about Gamma Ray Bursts. Computation of theoretical 

I Gamma Ray Bursts rate indicates that our models have too fast rotating cores, and that we need 

^ to include an additional effect to spin them down. Magnetic fields in stars act in this direction, 

and we show how they modify the evolution of massive star up to the final stages. 

-_. . Keywords, stars: evolution, stars: magnetic fields, stars: mass loss, gamma rays: bursts 

Q 1. Effects of rotation on the stellar surface 

-4— > Rotation has a strong influence on the steUar surface. Indeed, it adds a centrifugal 

—j ' component to the gravity, which modifies the shape of the surface, and various quantities 

1—^ , such as the effective temperat ure Tef f and the mass loss flux. These effects can be derived 
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from the von Zeipel theorem (jvon Ze ipel 1924), which is originally vahd for conservative 



cases of angular momentum dist ribution , and w as treated in the more general case of the 
so-called "shellular rotation" bv lMaedeiJ ( 19991 ). This theorem gives the relation between 



\^ , the local flux F of the star and the local effective gravity: 



where L{P) is the luminosity on the isobar and gofr the local effective gravity. The two 
remaining terms are given by 
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Here, p^ is the internal average density, M^, represents the effective mass, modified by 
the rotation velocity fi, x = AacT^ / {inp) is the thermal conductivity coefficient and xt 
is its partial derivative with respect to T, = ^ is the ratio of the horizontal density 



fluctuation to the average density on the isobar ( Zahnlll992f ). 5 is the thermodynamic 



coefficient 5 = — (9 In p/9 In T)p and Ht is the temperature scale height. Generally 
the term C,{0) is very small, and we can neglect it. 

The total gravity at the surface of the star is given by gtot = gcff + grad with giad = ^ 
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is the term due to radiative forces and k the total Rossland mean opacity. The Eddington 
limit in a rotating star is defined by the vanishing of gtot and we find the limiting flux 



-gcff . 



(1.4) 



The Eddington factor at a gi ven colatitude , which is given by the ratio of the local flux 
to the limiting flux, becomes ( Maederill999r ) 
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A first interesting consequence is that the maximum luminosity of a star, given by 
rn(^) = 1, is lowered by rotation (compare Lmax given above with the non-rotating 



one L 



max, no rot 



47TCGM \ 



Rot ation has also a strong impact on the mass loss rate of the star. lMaeder and Mevneti 
( 20001 ) found the following relation between the mass loss rate at a given rotational 
velocity M{n) and the non-rotating one: 
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where T — . '^j^,r is the "non-rotating" Eddington factor. We see that rotation will 
increase the total mass loss rate of the star. 

Rotation has another effect on the mass loss rate: it is no long isotropic at the su rface 
of a rotating star, but becomes colatitude-dependent. iMaeder and Mevneti pOOOT ) give 
the following local mass loss rate AM per unit surface Act: 
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with K and a the force multiplier parameters. We neglect here the small effect of C(^)- 
Rotation favors mass loss through the terms Mi, and Tq. Mass loss occurs preferentially 
where QeS is small, i.e at the poles. Then mass loss rate is thus varying as a function of 
the colatitude, producing the anisotropic wind phenomenon. 

Figure [1] shows various effects of the rotation at the surface of a 20 Mq star at a 
metallicity of 10^^ and at 95% of the critical rotation velocity. We can first remark that 
the star becomes oblate, with an equatorial-to-polar radius ratio -p-^3. ~ 13, Then we 
note the variation of the effective temperature with respect to colatitude: Toff at the pole 
is around 48000 K, while it is only around 34000 K at the equator. Finally, we see that 
the mass loss flux is larger at the pole by a factor of ~ 3.8. This allows the star to lose 
10% less angular momentum than if the same amount of mass was lost isotropically. 



2. Models without magnetic field and rate of GRB 

Let us briefly recall here the main assumpt ions of the so-called " collapsar" model for 
long-soft Gamma Ray Bursts tWo oslevlll993 ). In order to produce such an event, the 
following conditions must be fulfilled: 

• formation of a black hole; 

• enough angular momentum in the stellar core in order to form an accretion disk 
around the BH; 

• formation of a type Ic supernova (see lWooslev and BloomI ( 20061) ). 
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Figure 1. Effects of rotation on a 20 M© star at Z — 10 '' and with jf^ — 0.95. The star is 
seen equator-on, and the axis are in Rq. The latitudinal variation of Tcff is shown (gray scale), 
and the arrows represents the masse flux at a given colatitude. 

These three points agree with a fast rotating massive star. It is thus interesting to 
study the evolution of the rate of type Ic supernovae. Figure [2] shows the relative rates 
for type lb and type Ic SNe, computed with rotating models without magnetic field. 
To distinguish between type lb and type Ic, we use a criterion based on the amount of 
He ejected during the SN event: all models ejecting more than 0.55 M© of hehum and 
no hydrogen are considered as type lb, the models ejecting less than 0.55 Mq of helium 
(and still no hydrogen) are considered to give birth to a type Ic SN event. The interesting 
curve for our purpose is the dashed one, representing the SN Ic / SN II ratio with respect 
to the metallicity. Wee see that the number fraction of type Ic supernova becomes higher 
at higher metallicity. This is well in ag reement with the tr end shown by the observed 
type Ic to type II SNe ratios given by IPrieto et al.l (|2008l ) . There are many observed 
evidence s indicating th a t long -soft GRB occurs preferentially in metal-poor regions. For 
instance. iModiaz et al.l ( 20081 1 find that GRB events appear at low metallicity: between 
0.2 < Z/Zq < 0.7. Thus only type Ic events in metal poor regions (or part of them) can 
occur simultaneously with a GRB event. 

Moreover, if we consider our models at metallicities that are compatible with the 
observed GRB range of metallicity, we see that all the models producing a type Ic SN 
keep enough angula r momentum in the core to fulfilled the collapsar model conditions (see 



Hirschi et al. 



(J2005)). We can thus determine the ratio of GRB event to the total number 
of core collapse SNc; we fou nd that this rate is around 15% for metallicities between 0.4 

( 2004 ) found an observational rate of 



In comparison, iPodsiadlowski et al 



and 0.7 Zq. 

0.04% — 8%, depending of the aperture angle of the bipolar jets produced during the 
GRB event. Our theoretical rate is therefore much larger than the observational one. 
These two facts lead to the conclusion that not all type Ic SNe produce a GRB. We 
have thus to find a way to reduce the number of GRB progenitor candidates, in order to 
reproduce the observational rate. One possibility is to introduce new ingredients in our 
models to extract more angular momentum from the core during the stellar life. 



3. Models with magnetic field and wind anisotropy 

Following ISpruitI ( 20021 ). we have included in our models the effect of magnetic field 
amplified at the expense of the excess energy in the shear. As noted just above, this 
produces a strong coupling between the differentially rotating layers, and tend to build 
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Figure 2. SNIb/SNII (solid line), SNIb/SNII (dotted line) and SNIc/SNII (dashed line) ratios. 
The triangles are observed SNIb/SNII ratio at vari ous metallici t ies an d the upside down gray 
triangles observed SNIc/SNII ratio (extracted from lPrieto et alj (120081 ')'). 



a solid-body rotation profile. Contrarily to the previous models, where the rotational 
velocity is only weakly coupled between the surface and the core, models with magnetic 
field have a strong coupling, and thus, the loss of angular momentum due to mass loss is 
quickly transmitted to the core. This implies a strong extraction of angular momentum 
during the evolution, particularly when the mass loss is strong at the surface (e.g. during 
the Wolf-Rayet phases). 

To explore the combined effects of magnetic field and wind anisotropy, we computed 
two 60 M0 models a.t Z — 0.002 with fl/flci-it = 0. 75, with and without t he tre atment of 
the wind anisotropy. "With respect to the work by iMevnet and Maedeij ( 20071 ). we have 
improved significantly the treatnrent of the wind anisotropy, allowing this treatment to 
apply even when the critical velocity is reached and checking very carefully that the 
sum of the angular momentum remaining in the star and the angular momentum lost 
in the wind remains constant all over the evolution (see Fig. [3l bottom panel). As we 
shall see below, t his improvements lead to eff ects which although still important are less 
pronounced as in ' Meynet and MaedeiJ (J2007r ). 

For both models, the strong mixing induced by the magnetic field lea ds t o the so-called 

"quas i-chemically- homogenous" evolution fsee lYoon and Langeij (|2005l ) and lWoosley and Heger 
(^2006j). Figure[3] (top panel) shows the evolution of the total angular momentum of the 
stellar interior during its evolution. We see that the anisotropic model (dashed curve) 
is slightly higher than the isotropic one (solid curve). This is due to the effect of wind 
anisotropy, as shown in the medium panel: the anisotropic model lose less angular mo- 
mentum (around 7%), while the rotational velocity is high (first part of the evolution). 
Then, the star becomes a WR, and the high mass loss rate implies a strong breaking of 
the rotation: the surface velocity is sufficiently below the critical velocity for removing 
any anisotropy in the wind, and there is no more differences between these two models. 

We see that when magnetic field is accounted for, the WR phase has a crucial influence 
on the total angular momentum kept in the star (and thus in the core, through the 
coupling produced by the magnetic field). Our model does not keep enough angular 
monrentum for the coUapsar model, and is thus not a good GRB progenitor candidate. 
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Figure 3. Top panel: evolution of the total angular momentum kept in the star (in units of 
[10*^ gcTn^ s^^]). The dased curve is the model with anisotropic wind, the solid curve for the 
isotropic one. Middle panel: angular momentum removed by wind at each time step, in the 
same units: dased curve for anisotropic model, solid for isotropic one. Bottom panel: Sum of the 
angular momentum of the star and the integrated angular momentum removed by wind. 

However, models with very high initial rotational velocity would most probably develop 
larger differences between the iso- and anisotropic treatment. The same would be true 
for models with lower mass loss rate (less massive models, lower metallicity). We will 
adrcss this point in a forthcoming paper. 
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